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Abstract. We describe how the notion of optical beam shifts (including the spatial 
and angular Goos-Hanchen shift and Imbert-Federov shift) can be understood as a 
classical analogue of a quantum measurement of the polarization state of a paraxial 
beam by its transverse amplitude distribution. Under this scheme, complex quantum 
weak values are interpreted as spatial and angular shifts of polarized scalar components 
of the reflected beam. This connection leads us to predict an extra spatial shift for 
beams with a radially- varying phase dependence. 
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1. Introduction 

There are many analogies between phenomena in quantum theory and in classical wave 
optics, including the connection between Heisenberg's uncertainty principle and the 
bandwidth theorem p[], transverse optical polarization and spin half as 2-state systems 
[2] , and the Schodinger equation and paraxial equation of hght [3] , connecting the time 
propagation of quantum wavepackets and the propagation of narrow, coherent light 
beams. Although the underlying physics differs in these cases, the similarity in the 
underlying mathematics gives rise to analogous phenomena. 

Our aim in this paper is to describe and explore the strong analogy between 
the spatial and angular shifts that light beams experience on reflection from a planar 
interface, and the notion of weak quantum measurement, incorporating quantum weak 
values of operators [U El E] . Such a connection was identified by Hosten and Kwiat [7] 
in the experimental measurement of the spin Hall effect of refracted light beams via 
quantum weak values. Here, we generalize this approach, and show that the shift of 
a paraxial light beam, and its polarized components, is determined by the appropriate 
average value of an 'Artmann operator' related to the reflection matrix for the beam. 

Spatial beam shifts are usually small: they are typically comparable to an optical 
wavelength, which is much smaller than the physical width of the paraxial beam, even 
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when focused on the reflecting interface. Angular shifts of the direction of the beam, 
when they exist, are also small, of the order of the (narrow) spectral width of the beam in 
Fourier space. The effect originates physically from the Fresnel coefficients' dependence 
on incidence angle - as the plane wave components making up the beam deviate from 
their mean, the complex amplitude and possibly polarization of the reflected plane wave 
components varies. The magnitude and direction of the beam shift thus depends on the 
incidence angle of the beam and its incident polarization, and while we concentrate in 
this paper on reflection, the case of transmission is very similar. We will not derive 
all beam shift formulas here, as they involve some technical details unimportant to 
the connection to quantum measurements. For these details (in an explicitly optical 
approach), we refer the reader to a companion paper [8]. 

Historically, the first beam shift to be discovered was for light beams polarized 
linearly in, or perpendicular to, the plane of incidence {p- and s-polarizations 
respectively), by Goos and Hanchen in 1943 [9]. In this case, the shift is in the plane 
of incidence (referred to as 'longitudinal'), given by the famous Artmann formula fiU\ . 
If the incident beam has circular polarization (or, more generally, has any polarization 
not linear in or perpendicular to the plane of incidence), there is an 'Imbert-Federov 
shift' transverse to the plane of incidence, which has been explained in terms of spin- 
orbit coupling [TT]. Goos-Hanchen and Imbert-Federov shifts are largest in the regime 
of total reflection, close to the critical angle; when reflection is partial (reflection by a 
denser medium, or inside the critical angle), Fourier flltering of the spectrum leads to 
an angular shift of the propagation direction of the light beam fl2[ [T3t [H], which has 
been the subject of much recent experimental attention fiEl [T6] . 

A schematic showing the geometry of a longitudinal shift, both spatial and angular, 
is shown in Fig. [l| It is natural to introduce the concept of a 'virtual beam', represented 
in red in the flgure. Inspired by the method of images, this is an ideal, unphysical 
reflected beam whose position and direction is as if the reflection were purely specular. 
Rather than refer to spatial and angular shifts with respect to the incident beam, we 
instead denote the shifts with respect to the expected position and direction of the 
virtual beam. It therefore suffices to work in coordinates based on the virtual beam, 
which propagates in the ^-direction, and the plane of incidence is the xz-pleme, as shown 
in the figure. This avoids introducing complicated 3D geometry for the refiection with 
respect to the incident beam. 

Of course, the usual Goos-Hanchen and Imbert-Federov shifts apply to the 
centroid of the intensity of the light. Since the refiected beam also acquires a weakly 
inhomogenous polarization, analyzing the reflected beam with a polarizer reveals further 
shifts of polarized components, such as the transverse spin Hall effect of light [HI [7] 
of circular components when linear polarization is incident; in this case, each circular 
component undergoes a shift even when there is no net transverse shift to the overall 
beam. 

A particular interest of beam shifts is that the nature and magnitude of the resulting 
shifts are independent of the spatial proflle of the light beam provided the incident 
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Figure 1. Schematic illustration of a spatial shift Dr^ and angular shift A^, of a 
reflected beam in the plane of incidence. If reflection of the incident beam (green) was 
purely specular, one would observe the 'virtual beam' (red). However, due to angle- 
dependence of the reflection coefficients, the real beam is displaced both in position 
and propagation direction (blue). 

beam is axisymmetric and focused on the interface. When this condition is relaxed, for 
instance with a beam with an azimuth-dependent phase vortex exp(i£0) [l7l [18], there 
is an additional 'vortex-induced shift,' where the spatial shift involves a combination of 
the Artmann expressions for the spatial and angular shifts [I9l [16] . 

Each of the notions above corresponds to a similar notion in quantum measurement 
via a weak von Neumann interaction, and our purpose in this paper is to make these 
explicit. In the optical reflection case, the interaction is between the polarization and 
spatial structure of the beam, and the weakness follows from paraxiality (i.e. localization 
in Fourier space). Before describing the details of the analogy, we review the relevant 
notions from quantum measurements, largely following the description and notation of 
the very clear exposition of the quantum case by Jozsa [20] . 

2. Weak quantum measurement 

The weak quantum measurement of an observable A of a quantum system in initial 
state \ipi) is achieved by considering the entire system in a product state of the system 
to be measured, together with the state of the measurement pointer \ip), i.e. The 
interaction between the measured system and the pointer is weak, in the sense that the 
interaction hamiltonian ifjnt = gAp, acting at a single instant in time, where g <^ 1 is 
a small coupling constant, and p is the momentum operator generating translations of 
the pointer state \ip). The unitary evolution operator corresponding to the interaction. 
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in units with = 1, is 

exp{-iHi^t) ^ 1 - igAp, (2.1) 

following from the fact that the coupling constant is small, i.e. the interaction is weak. 

Now, the measurement may be made in the usual von Neumann sense, where the 
measured system's freedoms are 'traced over', by taking the (sub-) inner product with 
\ipi), giving rise to the expectation value 

Alternatively, the projection of the measured system in a certain postselected final state 
\ipi) may be considered, giving rise to the so-called 'weak value' of the operator [U [5] 

A. ^ (2.3) 

When A is hermitian, (A) is necessarily real, although it is not for more general A. 
However, for any A, the weak value Ay^ is usually complex-valued, and depends on both 
the pre- and postselected states. 

Furthermore, the weak value may take on very large values ('superweak' [2T]), 



particularly when \^p^) and are almost orthogonal (so the denominator of Eq. (2.3) 
becomes vanishingly small). We will denote the average value of A, whether the 
expectation or weak value, by the generally complex 

_ (mm 



(2.4) 



where {ip\ is {ipi\ or 

The final pointer state, after the weak interaction and possibly postselection, is 
therefore 

(V'|exp(-iifint)|^i)|<^) ^ (^l^/'i) (1 - i^ap) Iv?) ^ {ilj\tlJi) exp{-igap)\ip) . (2.5) 

The result of the measurement is a shift in the mean position of the pointer wavefunction 
{q\^p) = ^p{q) = \(f{q)\&^^'^\ in the position representation, where x is the wavefunction's 
phase. When a is real, this simply corresponds to a translation of (p{q) by the small 
amount ga; when a is complex, the expectation value of q is shifted by 

D = gRe{a) - (?Im(a) j dggV'(g) = gRe{a) + 2glm{a) {qx' (q)) , (2.6) 

where the second term is nonzero when the pointer wavefunction ip{q) at the moment 
of interaction has a varying phase (so the probability current J{q) = \f\'^x' is nonzero), 
and we have assumed the unshifted expectation of position (g) = 0. The second line 
here follows from integration by parts, assuming x is sufficiently well behaved; in [20] , 
the term proportional to glm{b) was found to be mdt{(f) by the continuity equation 
for probability and Schrodinger's equation. In Fourier space, the pointer's conjugate 
wavefunction ^(p) also undergoes a shift to its mean position [22], [6], by 

A = (7(p2)lm(a), (2.7) 
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proportional to the imaginary part of a and the width (variance) of the Fourier transform 
(p{p). We note that it is common in the quantum mechanical literature to refer to 
the weak interaction between measured and pointer systems as a 'weak measurement', 
referring to the small magnitude of the coupling constant g; the measurement is weak 
regardless of whether the weak or expectation value of the operator is being measured. 

3. The analogy between beam shifts and quantum measurement 

We are now in a position to draw the analogy between beam shifts and quantum 
weak measurements. As previously in classical optics analogies with weak values 
[221 E], the polarization of the beam, represented by the constant transverse vector 
E, is identified with the measured system, and the complex amplitude of the beam 
corresponds to the pointer. The measured system is therefore effectively a 2-state 
system (i.e. a 2-dimensional Jones vector), and the pointer wavefunction (a normalized, 
position-dependent complex amplitude) ip{r) depends on r = {x,y) in the plane 
perpendicular transverse to the propagation in z. A homogeneously polarized light beam, 
in the quantum language, therefore corresponds to the product state. Free paraxial 
propagation is analogous to free quantum time evolution according to the Schrodinger 
equation. The paraxial approximation applied to the beam follows from its strong 
localization around a mean propagation direction in direction (Fourier) space. A regular 
or weak measurement of the beam is made depending on whether the centroid of the 
total intensity of the final reflected beam, or a polarized component of it, are considered. 

When a plane wave encounters a planar dielectric interface, it is reflected [21]: 
it changes direction according to the law of specular reflection, and its polarization E 
undergoes reflection according to the appropriate Fresnel coefficients and Vp (evaluated 
at the incidence angle): the p-direction perpendicular to the plane of incidence, and 
the s-direction parallel to it, are distinguished as eigenpolarizations of the incident 
polarization E. The resulting polarization is therefore R ■ E, where R is a reflection 
matrix, diagonal in the s,p basis with entries given by the appropriate reflection 
coefficients. It is therefore natural to choose the beam coordinates to correspond to 
these eigenpolarization directions: x corresponding to p, in the plane of reflection, and 
y perpendicular to it (s polarization), as in Fig. [l} 

When the homogeneously polarized beam - rather than simply a plane wave - 
encounters a planar dielectric interface, it approximately undergoes the same change: it 
is specularly reflected, and its polarization changes according to its central wavevector 
component. However, the small variation of wavevector directions about the propagation 
direction gives rise to a small variation in the incidence angle and plane of incidence of 
each Fourier component, and it is this small variation, ignored in the approximation of 
beam as plane wave, that corresponds to weak quantum interaction. 

This is made explicit by the introduction of the virtual beam, that is, the imaginary 
beam propagating in the z-direction (i.e. after specular reflection) as discussed above. 
The virtual beam is simply specularly reflected, and its polarization is constant; the 
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'weak interaction' is then the residual small correction (to first order) from the angle 
dependence. It is convenient to perform the analysis in the coordinates of the virtual 
beam. 

Spherical angles for beam coordinates will be denoted by azimuth a and colatitude 
S, i.e. the transverse wavevector (whose length k is fixed) 

K = {Kx, Ky) = ksm6 (cos a, sin a) ^ k6 (cos a, sin a), (3.1) 

since we assume the spread of the beam in Fourier space is small: paraxiality implies that 
the variance (5^) = p- J cPK \K\'^\ip{K)\'^ <^ 1. The polarization of the virtual beam 
is assumed to be uniform, and determined by R, the refiection matrix of the mean 
wavevector. R is the refiection matrix for a plane wave incident at ^o, the incidence 
angle of the centre of the beam. 

The virtual beam is therefore in a product state, whose amplitude distribution is the 
same as the initial beam, and whose homogeneous polarization is given hj Ei = H ■ E. 
It is this virtual beam which plays the role of the preselected quantum product state: 
the main change to the beam on refiection is its specular change in direction and plane 
wave refiection-like change to its polarization: the smaller shifts in mean position and 
direction follow from the weak 1^-dependent variation in R. 

The physically refiected beam, in its Fourier representation, is found by 
multiplication of the virtual beam Fourier transform with the appropriate 1^-dependent 
refiection matrix R. Each K has its own particular plane of incidence, and R applies 
the appropriate refiection coefficients in the local s,p basis. Since the spread of K- 
directions around the propagation direction is small, the beam is sensitive only to a 
low-order Taylor expansion of this matrix R with respect to 6, that is 

R ^ R + 5( cos tt, sin a) ■ (R^,Rj,), (3-2) 

where R^; and R^^ represent the mean derivatives of R in the longitudinal and transverse 
directions respectively, evaluated at the mean wavevector. In this paper, we will not 
perform the extra geometrical calculations required to find explicit forms of R^ and R^; 
we derive these in the companion paper [8] in detail. 

It should be clear by analogy with the exposition in Section |2] above that in this 
approximation, the refiection operator may be rewritten directly as the action of the 
mean R followed by an interaction operator entangling the position and polarization. 
This can be represented by the mean refiection R left-multiplied by a weak interaction 
operator (as described above in Section |2]), 

R ^ (1 - iK ■ (A,, A,) ) R ^ exp {-iK ■ (A„ A,) ) R, (3.3) 

where (A^., A^) is a vector of 2 x 2 matrices we call 'Artmann operators', 

Aj = ^B~R'\ j=x,y. (3.4) 

If the refiection is total (i.e. the refiection coefficients are unimodular), it is 
straightforward to see (but omitted here) that each Aj matrix is hermitian. However, 
if the refiection is partial, some of the incident light is lost through transmission. 
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resulting in the Aj matrices being nonhermitian. The corresponding evolution operator 
is nonunitary as it does not preserve normalization, as some light is refracted. 

In Fourier space, the reflection operator acts as an impulsive evolution of the virtual 
beam under the 'interaction hamiltonian' K-{Ax, Ay), which is weak since the transverse 
momentum K is necessarily small for the paraxial beam. Apart from the generalization 
to two dimensions (which is straightforward unless (p{r) is complex, for which see Section 
|4]), the connection with the quantum measurement described above is immediate: the 
centre of the physical beam in real space, and its direction in Fourier space, are effectively 
a weak measurement-like shift to the pointer if{r), and its Fourier transform (f{K) by 
the real and imaginary parts of the (possibly complex-valued) average of the Artmann 
operators A^^, Ay. 

In the case of optical polarization, 'tracing out' the polarization degrees of freedom 
is simply considering the overall intensity of the beam ignoring the polarization; in this 
case, the spatial shift in j = x,y is given by the general shift (for a real <f{r)) 

-u 



^ Ei-Ei yk Ei-E; 
= -ImU: ^ . (3.5) 

When E is linearly polarized in the x- or y-directions (p and s polarized respectively), 
the corresponding longitudinal Goos-Hanchen shift is given by the famous Artmann 
formulas — ^Imr^/r^, — ^Imr^/r^ [10], which follow directly from the spatial shift formula 
in Eq. (3.5), and Dy = in this case. When E is circularly polarized (right- (+) or left- 
(— ) handed), there is a transverse Imbert-Federov shift Dy given by i^l^^s + rppcot^o 
[25| l26| fTTj [271 E]. It can be shown that Ry, originating from spin-orbit coupling, is a 
constant times the Pauli matrix ax in the s,p basis, so in fact Dy is zero if E is linearly 
polarized. 

If reflection is total so the Aj are hermitian, these spatial shifts are the only shifts 
which occur. However, if reflection is partial, there is a shift to the Fourier transform 
ip{K), which is k times an angular shift to the direction of the beam in the x and 
y-directions. From above, this is proportional to the variance of the beam's Fourier 
distribution (l-K'P) = k'^{6'^), where (S"^) is the variation of the beam in direction space. 
Thus 

^ ^ (W!>i,..^.--A,.iS, 
^ k Ei-Ei 

= {S')Re ( ^ ^5 ^ 1 (3-6) 
V E ■ R^R ■ E J 

This equation is derived purely from the optical viewpoint in [8]. Despite the recent 
interest in the angular shift, the simple form here in terms of the variance in direction 
space has not previously been emphasized. The universality, coming from the relation 
with weak interactions (although the value of the nonhermitian Artmann operator is 
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not weak here) is one of the main new statements in optical beam shifts in this paper. 

What of weak values themselves? As discussed above, postselection is clearly here 
represented by the presence of a polarizing analyzer F, which projects the beam onto a 
'final' polarization state. In terms of optical physics, the analogue of the entanglement 
between measured system and pointer is a position-dependent polarization pattern: the 
mean position of different final polarization states varies (Fig. [2]), depending on the 
details of the incident polarization and reflection matrix. The component shift, or 
'weak' shift, can therefore be written as follows: the component spatial shift is 

D^. = Re — — - = --Im =d , 3.7 

^ F-Ei k \F* -R-E J ^ ^ 

again for j = x, y, and the component angular shift is 

(IKR F*-A.-Ei fF*-B~-E\ , , 

^. = ^^im— — ^— ^ = 52 Re _J . 3.8 

^ k F* -Ei \F* - R-E J 

Of course, the weak shift typically has an angular contribution even when the are 
hermitian, representing the inhomogeneity of the polarization pattern in direction space, 
and weak shifts can be rather different directions to their strong counterparts (or exist 
when these are zero, as is the case of the spin Hall effect for light f7]). 

An illustration of optical beam shifts, including component shifts, is shown in Fig. [2] 
In each case, on the lengthscale of k~^, the inhomogeneous polarization is extremely 
weak, although both longitudinal and transverse shifts are present. A systematic 
theoretical study of the component and beam shifts for different choices of E and F is 
made in [8]. 

For certain choices of F, the component shift may be very large, with a rather 
small overall amplitude. However, it is important to note that such superweak shifts 
occur when F* ■ E^ = F* ■ R ■ E is small, and the inner product between initial and 
final polarizations F* ■ E does not directly play a role. This is of course because the 
analogue of the prepared quantum state is the virtual beam. 

These optical shifts completely mirror quantum mechanical pointer shifts, and the 
2-dimensional nature does not change any of the essential description. However, if the 
complex pointer amplitude (p{r) has radial- and azimuthal-varying complex phase, the 
situation is more complicated than the one-dimensional case described in [20], and we 
describe this in the following section. 

4. The case of radial and azimuthally varying pointer phase 

In this section, we will represent the weak or expectation value of the Artmann operator 
by the complex 2-dimensional vector a = {ax,ay), which is defined by one of the 
following: 

E • R Rj ■ E F* ■ Rj ■ E r a \ 

ai = =t= or p 3=x,y. (4.1) 

E R^R E F R E 
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Figure 2. Examples of shifted, reflected beams with weakly varying polarization, 
in units of k~^. In each case, the incident beam has polarization E = '^i^'^)^ with 
an incident gaussian profile with width lOOfc"^ at incidence angle = arcsin3/4. 
Reflection is between air and glass, with (a) total reflection (n = 3/2), (b) partial 
reflection (n = 2/3). The plots, in virtual beam coordinates, show contours of overall 
intensity (grey), and (a) magnitude of components in p- and s-polarization (purple 
and green respectively); (b) magnitude of right-hand and left-hand circular components 
(cyan and turquoise respectively). Despite the beam shifts being different, the variation 
in polarization from the virtual beam polarization Eq is imperceptible at this scale 
(blue ellipses and lines). 



Although our main example is optical beams, the discussion applies to any weakly 
interacting quantum system with two noncommuting operators Ax,Ay multiplying the 
two components of pointer momentum, and our aim is to generalize the contribution of 
Im(a) to the spatial pointer shift. 

Following [20], we consider the modulus squared (intensity or probability 
distribution) of the reflected wave after tracing out or postselecting the polarization 
state, 

I (1 - iK ■ a)(^(r) |2 ^ le-^^'^y^^l^ - i Im(a) ■ (^* V<^ - ^V<^*) (4.2) 

where K, as the momentum, corresponds to the generator of translations — iV, and only 
terms up to the first order in a are kept as usual. The first term represents the intensity 
pattern shifted by Re(a), and the second term represents additional interference caused 
between the differently-weighted Fourier components with varying incident phases. 

This second term generalizes the second term on the RHS of (2.6) to 2-dimensional 
pointer wavefunctions; the first term here - a pure translation - is just the usual spatial 
shift. The shift is the mean of r with respect to this distribution, and we assume that 
the mean position of the unshifted (virtual) beam is the origin. 
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We begin our analysis of the second term of Eq. (4.2) by rewriting it as the (optical 
or probability) current vector J = —^{(p*'Vip — (pVip*), 

- i Im(a) ■ {ip*Vip - (fVifi*) = 2 Im(a) ■ J. (4.3) 

Its contribution to the x-component of the spatial shift is therefore 

^extra ^2 J d'^r X Im(a) ■J = -J dV x^d^ [Im(a) • J] , (4.4) 

with the equality following from Green's theorem. An equivalent expression holds for 
the shift in the y-direction. 

In the one-dimensional case [20] , this is related through conservation of current and 
the Schrodinger equation to the rate of change of the variance of the wavefunction as in 
Eq. (2.6). This is not possible in two or higher dimensions since the integral does not 
involve the divergence of J, and the general situation appears rather more complicated. 

However, in most optical fields of interest, any varying phase factor factorizes into 
an azimuthal part (say with a vortex with quantum number £) and a radial part, i.e. 
for azimuthal angle and radius r. 



for intensity /(r) 



/rWexp(i£0 + i/(r)), 



(4.5) 



^, and /(r) some radius-dependent phase function. It should 
be noted that, for optical reflection, the azimuthal index i reverses sign on reflection 
[28|l29]. Our convention here is that Eq. ( |4.5 ) refers to the virtual beam, so the azimuthal 
dependence of the incident beam is exp(— ii 



It is easy to see that, for the virtual beam (4.5) on returning to cartesian cordinates. 



(4.6) 



Using this form in the formula for the extra shift in Eq. (4.4) in x, and after integrating 
out the azimuth 6, we have 



D 



extra 



-Im(a) 



drr (r/(r)[f (r) + 3r/"(r)] +3rV(r)/'(r),£[3r/'(r) -2/(r)]) 
= Im(a)-((r/'(r)),£), (4.7) 

where the second line follows from normalization and integration by parts, assuming 
/(r) is reasonably well-behaved. An identical argument for y gives 

D-*- = Im(a)-(-£,(r/'(r))), (4.8) 



so the net total shift in two dimensions, analogous to (2.6 ) for a pointer with the assumed 
form, is 

D = Re(a) + ((r/'(r)) + i^cra) ■ Im(a), (4.9) 

where (T2 is the second Pauli matrix. 

This is the form of the the spatial shift for complex a and complex (p{r). The 
£-dependent parts correspond to the 'vortex-induced shift' |28l |30l UHl [16]: the £-fold 
twisting of the phase fronts results in Im(ay), usually determining the angular shift in y, 
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contributing to the real shift in x weighted by i and similarly (in a manner preserving the 
sense of circulation), —ilm{ax) contributes to the shift in y. Apart from the azimuthal 
quantum number i, no other features of the incident beam contribute to this part of the 
shift. 

The other part, given by Im(a)(r/'(r)), directly generalizes the corresponding term 
in (2.6). In optical beams, a radially- varying phase is interpreted as curved phasefront 
(rather than a beam at its focus, where the phasefronts are flat and /(r) = 0), due, 
for instance, to the Gouy phase [31l |32]. This radial phase variation leads to a term 
proportional to the expectation value of phasefront gradient (r/'(r)) times the usual 
angular shift. 

A simple example to consider is a gaussian light beam (gaussian wavepacket) of 
width w propagating according to the paraxial equation, {p{r,z) = exp(— r^/[w^(l + 
iz/zR)])/(l + iz/zji), where 2;r = kw'^/2, the Rayleigh distance. In this case, {rf'{r)) = 
z/zji, implying the radial part of the extra shift is small when the beam is focused close to 
the interface; this result naturally is similar to the one-dimensional [20] constant related 
to the rate of expansion of the transverse waveform. To our knowledge, this additional 
'defocusing- induced shift' has not been experimentally measured, but is present as an 
additional feature of a more complicated amplitude pattern, corresponding to a more 
complicated quantum mechanical pointer wavefunction [331 1^ . In particular it is 
different to a focal shift, that is a shift of the focus along the propagation axis of 
the beam as discussed in 



5. Discussion 



We have shown how the reflection of a narrow optical beam is a classical wave analogue to 
a quantum weak measurement. Although the experiment could be done for an ensemble 
of single photons, nothing in our analysis relies strictly on the quantum nature of the 
light: the Hilbert space is completely classical, describing polarization and (transverse) 
position of the optical beam, with entanglement of the degrees of freedom interpreted 
as a position-dependent polarization pattern [55] . 

As such, this phenomenon adds to many well-known classical optical analogues of 
quantum phenomena involving polarization and complex optical amplitudes. Although 
our paper highlights the analogy between the weak values and optical beam shifts of 
polarization components, the discussion applies to the whole beam without a polarizing 
analyzer, which is itself analogous to a usual expectation value. In particular, the 
effect of the weak measurement on the imaginary part of an operator's average value 
is immaterial of the kind of average (weak or expectation); the shift to the momentum 
wavefunction and possibly extra spatial shift always occurs for average values of 
nonhermitian operators, realised here in partial reflection. 

The main effect of reflection, of course, is to change a beam's propagation direction 
and polarization, giving what we call the virtual beam. It is the virtual beam which 
is analogous to the quantum prepared system, whose position and polarization weakly 
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interact, causing it to be shifted in both position and direction. Rather than coming 
from a small coupling constant, the weakness comes from paraxiality, that is narrowness 
of the beam about its mean propagation direction in Fourier space. In this case, as we 
have seen, the main change to a beam on reflection is the small shift in position and 
possibly direction; for nonparaxial beams (such as the dipole radiation field considered 
in [36]), the shifts are present, but the reflected field is much more complicated. 

The effects described here are only to first order, as are most phenomena studied in 
the physics of weak measurement. Second order effects become important for superweak 
values, whenever the analyzer is almost orthogonal to the polarization of the virtual 
beam. The magnitude of the spatial and angular shifts corresponding to superweak 
values is explored in an optical setting elsewhere [S^- With the exception of the 
extra shift coming from the pointer wavef unction's varying phase, all of the shifts are 
independent of the spatial distribution of the amplitude of the beam (assuming its 
modulus is radially symmetric). This simplifies previous descriptions of the angular 
shift in optics, as the shift has a universal form when measured in units of the angular 
spread (5^). Higher-order effects involve changes to the shape of the beam, and these 
will reveal more subtle structures associated with the entanglement between position 
and polarization. 

It should be stressed that the shift of a beam on reflection does not require 
polarization, only an incidence angle-dependent reflection coefficient r{6o). This might 
be achieved with acoustic waves at a lossy wall (similar to the partial reflection case), or 
Robin boundary conditions (similar to total reflection) |38]. From the discussion above, 
it is clear that even in this case there is a complex scalar shift a = (i//c)(r'/r), whose 
real and imaginary parts contribute in the ways described above: the real part gives rise 
to a longitudinal spatial shift, the imaginary part to a longitudinal angular shift, and 
possibly to an additional longitudinal spatial shift (if the wave has a radius-dependent 
phase) and a transverse shift (if there is an azimuthal phase). Most of the transverse 
shifts described above are absent as they require spin (polarization) interacting with 
azimuthal (orbital) terms. 

The fact that weak values generalise so simply to classical wave physics cements 
their significance as physical quantities. Furthermore, the generality of the beam shift 
framework as we have outlined suggests generalization to other kinds of waves, not only 
acoustic waves as suggested above, but also elastic waves and shifts to matter waves 
[391 1^ . This suggests that further analysis of such simple physical phenomena as the 
reflection of a beam at a dielectric might reveal new insights into quantum physics. 
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